Trapping of three-dimensional Holstein polarons by various impurities 
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We study the bound states of a three-dimensional Holstein polaron near various kinds of single 
impurities, using the momentum average approximation. We show that the electron-phonon coupling 
is responsible for a strong renormalization of the impurity potential, resulting in an effective potential 
with significant retardation effects, which describes essential physics ignored by "instantaneous" 
approximations. The accuracy of our approximation is gauged by comparison with results from 
Diagrammatic Monte Carlo for the case of an impurity that modifies the on-site energy of the 
electron. We also discuss impurities that modify the local strength of the electron-phonon coupling, 
as well as isotope substitutions that change both the electron-phonon coupling and the phonon 
frequency, and contrast and highlight the difference between these cases. 
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I. INTRODUCTION 

The challenge to understand the effects of disorder 
on the behavior of particles strongly coupled to bosons 
from their environment is commonly encountered in cor- 
related electron systems. For example, high-T c cuprates 
are doped anti-ferromagnetic insulators in which, be- 
side strong coupling to magnons, ARPES measurements 
have also provided evidence for strong electron-phonon 
coupling X At the same time, charge carriers moving in 
the Cu02 layers are subject to random disorder poten- 
tials from the adjacent dopant layers. Substituting only 
a few percent of the Cu atoms with impurities sup- 
presses superconductivity by localizing the low energy 
electronic states^ Inhomogeneities in the superconduct- 
ing gap measured in high resolution tunnelling exper- 
iments have been attributed to atomic scale disorder 
in the phonon energy and the electron-phonon coupling 
strength in these materials^ Organic semiconductors are 
another class of materials where interplay between dis- 
order and electron-phonon coupling is believed to be im- 
portant in determining their properties, and are currently 
under active investigation^^ 

Although the results we present here are valid for any 
type of bosons (so long as they can be modelled as dispcr- 
sionless Einstein modes), for simplicity in the following 
we restrict our discussion to optical phonons. The re- 
sult of the interplay between disorder and coupling to 
such phonons depends on their relative strengths. Dis- 
order that is considered weak for free electrons can be 
strong enough to localize a polaron, that is the dressed 
quasiparticlc which consists of the electron together with 
its cloud of phonons, because of its heavy effective mass. 
On the other hand, whereas in the weak disorder regime 
electron-phonon coupling hinders the motion of electrons, 
such coupling actually facilitates the electron mobility in 
the strongly disordered regime where the Anderson lo- 
calization prevails^ 

Certain aspects of this problem have been studied with 
various approximations, most of which rely on sophisti- 



cated non-perturbative methods^ such as the statistical 
dynamic mean field theory (DMFT)j&£ or dynamical co- 
herent potential approximation (DCPA) ^Siii The under- 
lying meaning of these approximations and their accu- 
racy is rather hard to gauge. On the computational side, 
refined versions of the approximation-free diagrammatic 
Monte Carlo (DMC) technique^ and of the continuous 
quantum Monte Carlo algorithm^ have recently been ap- 
plied to the problem of a Holstein polaron near a single 
impurity. While essentially exact, such calculations re- 
quire significant computational resources, and cannot be 
easily generalized to other couplings, for example. 

Here, we study the bound state formation for a three- 
dimensional Holstein polaron in the presence of an im- 
purity, using a generalization of the Momentum Average 
(MA) approximation to inhomogeneous systems J^ MA 
was originally developed to study homogeneous systems 
with various types of electron-phonon coupling J^r— It is 
a non-perturbative method that sums in a closed-form 
expression all the self-energy diagrams up to exponen- 
tially small corrections that are neglected. The method 
can be systematically improved,— therefore providing a 
fast yet accurate way to scan a vast range of parameters. 
The method can be used to study all possible types of 
disorder for various types of electron-phonon coupling. 
Here, we use it to consider different types of disorder 
due to single impurities, namely a variation in the on- 
site energy, in the electron-phonon coupling and/or in 
the phonon energy are separately considered. The accu- 
racy of this method is demonstrated for the former case 
by comparison with available DMC results. 

Unlike most other theoretical approximations, MA has 
the important benefit that its structure reveals the es- 
sential physics of such problems. It is well-known that 
electron-phonon coupling leads to the dressing of the par- 
ticle, resulting in a polaron with a larger effective mass. 
What MA reveals is that the electron-phonon coupling is 
also responsible for a renormalization of the disorder po- 
tential. This renormalization can be very large and has 
strong retardation effects. Moreover, the renormalizcd 



potential can have a finite-range even if the bare disor- 
der is on-site only. The single impurity problem provides 
us with a simple test case to understand the effects of this 
renormalization, and to accurately compare and contrast 
the behavior of the polaron in the presence of various 
types of local disorder. Such results are a necessary first 
step in order to gain the intuition needed for understand- 
ing the behavior of more complicated systems. 

The paper is organized as follows: In Section II we 
describe our method and discuss its meaning (full details 
are provided in the Appendix). Section III presents our 
results for the three types of impurities, and Section IV 
contains the summary. 



II. MOMENTUM AVERAGE APPROXIMATION 
FOR INHOMOGENEOUS SYSTEMS 

For completeness, we present the MA formalism for 
the general case of random on-site disorder plus inhomo- 
geneities in both the coupling and the phonon frequen- 
cies. A simpler case (with only on-site disorder) has been 
briefly discussed in Ref. [lj. The Hamiltonian is: 



n = n d + v el - ph =n Q + v d + v e i- ph 



(1) 



where T-Ld describes the non-interacting part of the 
Hamiltonian, and for convenience is further divided into: 

n Q = -tJ2(4 c j + h - c -) + J2 Qib t bi 

which contains the kinetic energy of the particle and the 
energy of the boson modes {% = 1), plus 



v<i = y^ tic\ci 



describing the on-site disorder. The interaction part: 

^„.„t„/At 



Vei- P h = 2_^ g t c}ci{b] + hi 



describes the (possibly inhomogeneous) Holstcin-like 
coupling^ between the particle and the bosons. Here, i 
indexes lattice sites - the lattice can be in any dimension, 
and of finite or infinite extent. The operators c,; and hi de- 
scribe, respectively, particle and boson annihilation from 
the corresponding state associated with lattice site i. The 
spin of the particle is ignored because it is irrelevant in 
this case, but generalizations are straightforward^ For 
simplicity, we assume nearest-neighbour hopping; this 
approximation can also be trivially relaxed. Depending 
on the model of interest, the on-site energies e^, electron- 
phonon couplings gi and phonon frequencies f2^ can be 
assumed to be random variables. As detailed below, our 
results here will focus on single on-site impurities such 
as €i = —USifi, but the formalism applies for any model 
consistent with Eq. ([TJ. 



Our goal is to calculate the single polaron Green's func- 
tion in real space: 



G iS (w) = (0\*&(w)ct\0) 



(0| Ci |a)(a| C t|0) 



^' uj-E a 



"I 



(2) 



where |0) is the vacuum, G(w) = [u — H + if]}^ 1 is 
the resolvent with r\ — > 0+, and E a , \a) arc single po- 
laron eigenenergies and eigenstates of the Hamiltonian, 
W\a) — E a \a). Knowledge of this Green's function al- 
lows us to find the spectrum from the poles, and the local 
density of states (LDOS) measured in scanning tunneling 
microscopy (STM) experiments, p(i;uj) = — —ImGaiu)). 
To calculate this quantity, we use repeatedly Dyson's 
identity: G(w) = G d (oj) + G(oj)V et ^ ph G d (uj), where 
G d {uj) = \u) — T-Ld + ir)} -1 is the resolvent for the non- 
interacting system. The first equation of motion gener- 
ated this way reads: 



G«( W ) = g£.(w) + S^M^M 



i 



where 



GUw) = (0| Ci G d ( W )4|0) 



(3) 



(4) 



arc, in principle, known quantities and we have intro- 
duced the generalized Green's functions: 



?{n) 



H = (0| Cl GHc]6]"|0). 



,(0), 



Note that F-j (w) = Gy(w). Next, we generate equa- 
tions of motion for these generalized Green's functions. 
For any n > 1, we find: 

F^M = 5>G&.( W - n%)(0|c,GM C t&^"|0) 

f^h + ^-^hI (5) 



+ i G>£.(w-nfi i ) 



This equation relates F^ n > not only to other Green's 
functions of a similar type, but also introduces new 
propagators with phonons at two different sites. Equa- 
tions of motions can be calculated for these new general- 
ized Green's functions, linking them to yet more general 
Green's functions, and so on and so forth. The resulting 
hierarchy of coupled equations describes the problem ex- 
actly, but is unmanageable. Approximations are needed 
to simplify it and find a closed-form solution. 

The main idea behind the MA approximations is to 
simplify this set of equations by neglecting exponentially 
small contributions in each equation of motion. At the 
simplest level - the so-called MA' ) approximation - we 
ignore the first term in Eq. $5§ for any n > 1. This 
is reasonable at low-energies like the ground state (GS), 
lu ~ Eg Si where u> — nQj is well below the energy spec- 
trum of Hd and, therefore, GfAu — nflj) is guaranteed 
to decrease exponentially with increasing distance \l — j\. 



As a result, here we keep the largest / = j propagator, 
and ignore exponentially smaller / ^ j contributions. Al- 
though this is the simplest possible such approximation, 
it is already accurate at low energies, as shown in the re- 
sults section. It can also be systematically improved, as 
discussed below. First, however, we complete this MA^ '- 
level solution. 

The simplified equation of motion now reads: 



F<?\u>) = 9j G^Aw - nil,) \f£ +1 \uj) + nF^ 1 ^) 



On physical grounds, we know that F^" (to) must vanish 
for sufficiently large n, because its Fourier transform is 
the amplitude of probability that a particle injected in 
the system will generate n phonons in time t, and this 
must vanish for large enough n. Asa result, these recur- 
sive equations admit the solution: 



4; i) h 



A n (j,u)F£- l) (uj), 



where the continued-fraction: 



A n (j,w) = 



ngjGjAu - nQ.j) 



l-gjGUu-nSl^An+i&u) 



(G) 



can be efficiently evaluated starting from a cutoff 
An c (j,uj) = for a sufficiently large N c . Generally 
speaking, this cutoff N c must be much larger than the 
average number of phonons expected at site j; in prac- 
tice, the cutoff is increased until convergence is reached 
to within the desired accuracy. Substituting F^ ■ (w) = 
A\(j, (jj)Gij(ui) in Eq. ^ leads to a closed system of lin- 
ear equations for the original Green's function: 



G«(w) - G?» + Y,Gu(u>)giA 1 (l, U )Gf j (u>). (7) 



This equation has a similar structure to the equation 
linking the disorder Green's function to the free parti- 
cle propagator G\j'(uj) = (0|ci[w + if] - Ho] _1 c]|0) (in 
the absence of coupling to phonons), which is depicted 
diagrammatically in Fig. QJa), and which reads: 



GUu>) 



<#>(") 



E 
i 



GiMziGfh") 



(8) 



This analogy shows that coupling to phonons renormal- 
izes the on-site disorder e; — > e; + giA\(l, u>). Note that 
Ai(l,u>) depends not only on the local phonon frequency 
fi; and coupling gi, but also on all the bare on-site en- 
ergies 6^ through the disorder propagators G d . This is 
the simplest example of the emergence of a renormalizcd 
potential for this problem, that is made very transparent 
within the MA approximation. 

While Eq. ([7]) can be solved directly for a finite-size 
system, we can improve its efficiency and reveal a dif- 
ferent physical interpretation by explicitly removing the 
"average" contribution due to the electron-phonon inter- 
actions. Let g and £1 be the average values of the gi, Qi 



(a) 



(b) 
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FIG. 1. (color online) Diagrammatic expansion for (a) the dis- 
order Green's function Gfj(u)) (bold red line); (b) the polaron 

Green's function in a clean system G\Mdj) (double thin line), 
and (c) the "instantaneous" approximation for the polaron 
Green's function in a disordered system, Gfj (u>) (double bold 
red line). The thin black lines depict free electron propaga- 
tors, the wriggly lines correspond to phonons, and scattering 
on the disorder potential is depicted the dashed lines ending 
with circles. See text for more details. 



distributions. We assume that the on-site energy aver- 
age e = (a finite value results in a trivial shift of all 
energies). Then, let: 



A n (uj) 



nggo(uj — ntt) 



1 - gg (u - nQ)A n+ i(co) 



(9) 



be the continued fractions corresponding to these average 
parameters, where we use the short-hand notation: 



ffoM = G|f( W ) = i^ 



1 



£k + IT] 



for the on-site free propagator (in the absence of disorder, 
this quantity becomes independent of the site) . It is given 
by the momentum average of the free propagator, where 
Ek is the free-particle dispersion. 

The "average" renormalization of the on-site energy 
is now recognized to represent the corresponding MA' ' 
self-energy for a "clean" system, i.e. a homogeneous sys- 
tem with average coupling and phonon frequency: 

s ma(">M =gAi(w) 

see, for instance, Eqs. (II) and (12) of Ref. [la . 
Introducing the effective disorder potential: 

vo(l,u) = 9iAi{l,u) - £ MA (o)(w), (10) 

Eq. ([7]) can be rewritten as: 

G i3 iw) = G% J (Q) + J2G«WMl,u)Gi i (fi), (11) 



where Cj — uj — E MA ( ) (uj). This energy rcnormalization, 
uj — > Cj, reflects the fact that processes describing the 
formation of the polaron in the "clean" system have been 
explicitly summed. 

Besides being numerically more efficient, since now 
Vo(j,ui) contains only the fluctuations from the (not nec- 
essarily small) average value included in Cj, Eq. (fTOj) re- 
veals a different interpretation for the effects of the in- 
terplay between disorder and electron-phonon coupling. 

Consider first the meaning of GfAui), which would be 
the solution if we could ignore Vo(j, ui). In the absence 
of on-site disorder this term equals G), (Cj), i.e. the ex- 
pected solution for a polaron in the clean system, de- 
picted diagrammatically in Fig. QJb) (of course, the exact 
self-energy is here approximated by £ MA (o) (w)). Com- 
paring Figs. QJa) and[ljb), it follows that GfAui) is the 
sum of the diagrams shown in Fig. [He). 

At first sight, this seems to be the full answer for 
this problem, since these diagrams sum the contribu- 
tions of all the processes in which the polaron scatters 
once, twice, etc., on the disorder potential. This is cer- 
tainly the answer obtained in the limit of an "instan- 
taneous" approximation valid when f2 — ► oo, i.e. when 
the ions are very light and respond instantaneously to 
the motion of electrons. In this case, one can perform a 
Lang-Firsov transformation and after an additional av- 
eraging over phonons, one obtains an approximative ef- 
fective Hamiltoniani 13 ' 21 



(a) 



7-Linst = -t* y^x c t c 3 + h - c -) +y^ ( e 
(hj) i 



<£_ 



C Ci 



(12) 



where t* = te « a is the renormalizcd polaron hop- 
ping, and —g 2 /Q is the polaron formation energy (for 
simplicity, here we assume that the phonon energies 
and electron-phonon coupling are homogeneous, gi —> 
g,Qi — > fi). The Green's function of this Hamiltonian 
is also given by Fig. [He), if the polaron propagator is 
approximated by: 



<#>(*)-» 



N ^ 

k 



= ik(Ri-Rj) 



El 

n 



uj 



where e£ is the renormalized kinetic energy. Of course, 
using the full expression of G\MCj) is preferable since the 
self-energy E MA ( > (uj) describes much more accurately 
the overall energy shift and effective mass renormaliza- 
tion than those asymptotic expressions, besides also in- 
cluding the quasiparticle weight. 

That Gfj (Cj) cannot be the full answer becomes obvious 
if we consider what happens when we rewrite the clean 
polaron propagators in terms of free particle and phonon 
lines, i.e. we substitute the expansion of Fig. [ljb) in 
[ljc). Doing so reveals that within this "instantaneous" 
approximation, scattering of the electron on the disorder 
potential is allowed only when no phonons are present, 
see zoom- in in Fig. [IJc). However, we know that for 



i i 

i t i i = + . 



(b) 



_, L 



* * * 



+ Vt W 'isS + h' ?#S •r-S + Efc W *. 



FIG. 2. (color online) (a) Diagrammatic expansion for the full 
MA solution Gij (ui) (thick dashed blue line) in terms of the 
clean polaron Green's function (double thin line) and scatter- 
ing on the renormalized disorder potential e*(u>), depicted by 
dashed lines ended with squares, (b) Diagrammatic expansion 
of t*(ui). For more details, see text. 



moderate and large electron-phonon coupling, the prob- 
ability to find no phonons in the system is exponentially 
small, therefore the processes summed in Fig. QJc) have 
very low probabilities. 

What is missing in Fig. QJc) are diagrams describing 
the scattering of the electron on the disorder potential 
in the presence of the phonons from the polaron cloud. 
Their contribution is included through the renormalized 
potential vq(1,uj) in the second term of Eq. (fTTj) . Indeed, 
the full MA solution shows that the polaron scatters not 
on the bare disorder ei , but on the renormalized disorder 
potential 



e*(w) = ei +v q (1,uj), 



(13) 



as depicted in Fig. [2ja). The diagrammatic expansion of 
the additional term v$(1,uj), shown in Fig. [U(b), verifies 
that it indeed describes the effective scattering in the 
presence of arbitrary numbers of phonons. 

Taken together, the diagrams summed in Fig. [2] repre- 
sent all possible contributions to the polaron propagator 
in the disordered system. The MA' ' approximation con- 
sists in discarding exponentially small contributions from 
each of these diagram, as already discussed. MA' ' also 
has an exact variational meaning, namely of assuming 
that the polaron cloud can have phonons only on a single 
site, in direct analogy with the clean casei 18 ' 22 It is quite 
remarkable that all the diagrams corresponding to this 
variational approximation can still be summed analyti- 
cally in closed form, even in the presence of disorder. 

As is the case for the clean system, MA can be sys- 
tematically improved by keeping more contributions to 
Eq. ([5]). In particular, at the MA' 1 ) level, we also treat 
the equation for the F^ 1 ' functions exactly, and make 
the MA approximation of discarding exponentially small 
off-diagonal propagators only for n > 2. The logic here 
is that the propagators appearing in F^ 1 ' have the high- 
est energy, therefore the slowest exponential decay. The 
MA' 1 ) equations can also be solved in closed form. The 
details are presented in the Appendix. The final solution 



looks identical to Eq. (fTTj) . except the renormalized en- 
ergy is now £) = lo — E MA (ij (oj) while the renormalized 
potential i>o(/,w) is replaced by a more complicated, yet 
more accurate expression vx(l,u)). The meaning of all 
these quantities, however, is the same. 

To summarize, MA reveals that the role of electron- 
phonon coupling is two-fold. On one hand, it renor- 
malizes the quasiparticlc properties due to polaron for- 
mation, just like in a clean system (as revealed by the 
explicit appearance of the "clean" system self-energy). 
However, this coupling also renormalizes the disorder po- 
tential experienced by the particle, ti — > e*(w). As we 
show next for various types of disorder, this renormaliza- 
tion is non-trivial in that it has strong retardation effects, 
and has significant consequences. "Instantaneous" ap- 
proximations completely ignore this rcnormalization, and 
therefore miss essential physics. To be fair, in practice 
the "instantaneous" approximations are usually imple- 
mented in an improved, variational for m 21 i 23 which is cer- 
tain to be much more accurate. However, this leads to the 
necessity to calculate the variational parameters through 
a self-consistent loop, making the improved version com- 
putationally as complicated as DMFT and DCPA, which 
also have self-consistency loops. 

In contrast, MA gives a closed analytical expression for 
all quantities of interest, in a formulation whose meaning 
is very transparent, and whose accuracy can be system- 
atically improved. 



The LDOS at the impurity site is then found to be: 



p(P;u) 



1 



ImGU") 



Po(0;w) 



\l + UG^(w) 



where po(0;u>) = — -ImG 00 (ui) is the LDOS in the clean 
system (equal to the DOS, because of translational invari- 
ance). As a result, a bound state below the continuum, 
signalled by a delta- function peak in p(0; lo), occurs if and 
only if the denominator of Eq. (fT6| vanishes. For a 3D 
simple cubic lattice this means that an impurity bound 
state appears if there is an energy E < —Qt such that 
ReGg (E) = —1/U (below the continuum the imaginary 
part of G 00 (E) vanishes). This equation can be solved 
graphically to find that a bound state appears for any 
U * U ° = -ReT^M^ - 3.96*. 

In the presence of electron-phonon coupling the equa- 
tions are more complicated, but the idea is the same: we 
calculate the LDOS at the impurity site and compare it 
to the DOS of the clean system. If the former has a peak 
below the threshold of the latter, then a bound state ex- 
ists at that energy. We then vary U to find the critical 
value above which a bound state is guaranteed. More 
details about the impurity state, such as its localization 
length, statistics for the phonon cloud, etc., can be ex- 
tracted from the LDOS at sites in the neighbourhood of 
the impurity. Here we focus on identifying when bound 
impurity states are stable. 



III. POLARON NEAR A SINGLE IMPURITY 

We now apply the general formalism described above 
to the simplest type of "disorder", namely an otherwise 
clean 3D simple cubic lattice with a single impurity. Im- 
purities which modulate the on-site energy, the strength 
of the electron-phonon coupling and/or the phonon fre- 
quency are separately considered. We investigate under 
what conditions such impurities can trap the polaron. 

As a reference case, we first review here briefly the 
solution in the absence of electron-phonon coupling. In 
this case, the impurity can only modulate the on-site 
potential, e^ = —USi^, and the Hamiltonian reduces to: 



H d = -t 



<ij> 



(c}cj + h.c.) - Uclco =n Q + V d . (14) 



For this form of Vd, Eq. ([5]) reads: 

G?» = Gg>M - UGJ (u)G$(u>), (15) 



and is trivially solved to find: 



GiM = Gfhu)-U 



^)t,,\n(v)< 



1 + UG ° \oj) ' 



(16) 



Of course, because of translational and time reversal 
symmetry, Gg>(w) = G^. (w) = Gf\ (uj), etc. 



A. Impurity changing the on-site energy 

The Hamiltonian describing this case is: 

H =T + flJ2blh + gJ2 c Ui(h + bl) - Ucl C(h (17) 



where T is the electron's tight-binding Hamiltonian. We 
are interested in the attractive impurities, U > 0, when 
an impurity state can be bound near the impurity site. 
To find the LDOS at the impurity site, we need to solve 
Eq. (fTTj) to find Goo(w). Note that now Gfj(u) is known, 
being given by Eq. (fTB")). The free- particle propagators 

(0)/ N 1 ^ „ ik ( R i 



V ^ 



-, where for the simple cubic 



2tJ2i=i cosfci, can be calculated by doing 



lattice £k 

the integrals over the Brillouin zone. A more efficient 

approach, which we use, is discussed in Ref. [24|. 

In this case, the renormalized impurity potential de- 
cays fast at sites away from the impurity, because 
Gf^to) -> G z ( j 0) (w) when |2| -> oo, so that A x (l,u)) -> 
Ai{uj). Physically, this is because the impurity has less 
and less influence at sites far from where it is located. 
Mathematically, this follows from Eq. (|16p and the fact 

that G l0 (lo) decreases exponentially with the distance 
between site I and the origin, at energies below the free 
particle continuum, of interest here. As a result, in Eq. 




7.548 
co/t 



FIG. 3. (color online) LDOS at the impurity site p(0, u>) in 
units of 1/i, vs. the energy co/t, for (a) MA' ' with l c — 
for U = 1.8, 1.9, 1.95 and 2.0. The dashed line shows the 
DOS for the clean system, times 100; (b) MA (0) at U/t - 2, 
and cutoffs in the renormalized potential l c = 0, 1, 2. Panels 
(c) and (d) are the same as (a) and (b), respectively, but 
using MA (1) . Panel (e) shows DMC results from Ref. 0, 
for same parameters as (a) and (b). Other parameters are 
Q = 2t, A = 0.8, r)/t = 10" 3 . 



()11|) we only need to sum over sites I close to the impurity, 
and the system can be solved very efficiently. 

The appropriate value for this cutoff varies depending 
on the various parameters, but generically it decreases as 
the energies of interest become lower. Again, mathemat- 
ically this is due to the exponential decrease of the free 
particle propagator with distance, and the fact that this 
decrease becomes faster as w — > — oo. Physically, this can 
be understood as follows. First, let us explain why is the 
renormalized disorder potential non-local, even though 
the bare impurity potential is local. The answer is pro- 
vided by the diagrams which contribute to it, see Fig. 
H(b). Consider, for simplicity, the MA 1 - ) approximation. 
In this case, all phonon lines appearing in these diagrams 
start and end at the site I for which Vq{1, lo) is being calcu- 
lated - this is the site where the phonon cloud is located. 
However, the electron is found with various probabilities 
away from the polaron cloud, so it can scatter on the im- 
purity if this is located within the "radius" of the polaron, 
where the electron resides. In other words, the range of 
the renormalized potential is controlled by the polaron 




" a = .5t (MA ) 

— ■ Q. = 2t (MA°) 

£2 = 4t (MA°) 

— • « = 8t (MA°) 

v n = .5t(DMC) 

□ n = 2t(DMC) 

n = 4t (DMC) 

O n = 8t(DMC) 



FIG. 4. (color online) Phase diagram separating the regime 
where the polaron is mobile (below the line) and trapped 
(above the line). The effective coupling is A = g 2 /(6tQ), and 
the critical trapping potential U c is shown for several values 
of Q/t. The MA results (filled symbols) compare well with 
the DMC results of Ref. [J2 (empty symbols). 



size. From medium to large couplings, as the polaron be- 
comes smaller, the renormalized potential becomes more 
local. At small couplings, though, the distance between 
the electron and its phonon cloud can be appreciable, 
and the range of the potential increases. 

In Fig. 03a) we plot p(0,w) over a wider energy range, 
for several values of [/, within MA 1 - ). The dashed line 
shows the DOS of the clean system, multiplied by 100 
for visibility. For U/t = 1.8, 1.9, the impurity attraction 
is not sufficient to bind a state below the continuum, al- 
though the LDOS is pushed towards the lower band-edge. 
For U = 1.95£, there is a peak just below the continuum. 
Because of the finite value of -q, the two features are not 
completely separated and the continuum onset looks like 
a "shoulder" , however lowering r\ allows us to clearly sep- 
arate the two features (not shown). Finally, for U = 2t 
the bound state peak is clearly below the continuum, so 
in this case U c « 1.95£. This U c value equals that ob- 
tained in DMC, although our energies are overall higher 
than the exact DMC values shown in Fig. [3(c), as ex- 
pected for a variational approximation. At the MA' 1 ) 
level the agreement with DMC is significantly improved 
since all features move toward lower energies, see Fig. 
[3(c). The critical value U c ~ 1.95t at which an impu- 
rity state appears below the continuum is little affected, 
however, by this overall shift of the spectral weights. 

The dependence of the bound state energy on the cut- 
off is shown in Figs. [3(b) and (d) for MA'°) and MA' 1 ), 
respectively, for the case with U = 2t, O = 2t and 

2 

the effective coupling A = ^ = 0.8. At these ener- 
gies, keeping only the local part in MA' ), i.e. setting 
vq(1,lo) — > 5iflVo(l,u)), is already a very good approxima- 
tion. Including the correction from the 6 nearest neighbor 
sites (l c = 1) lowers the energy somewhat, but the con- 
tribution from the second nearest neighbor sites (l c = 2) 




FIG. 5. (color online) Effective value of the impurity attrac- 
tion U*/U, extracted from the scaling E B /t* = f(U*/t*), for 
A = 0.5, 1.5 and fi = 3t. In order to have the best fit to the 
data, we plot each curve starting from slightly larger U/t than 
the corresponding U c /t given in Fig. Q. For more details, 
see text. 



is no longer visible on this scale for either the MA' ' or 
the MAW results. 

Repeating this process for other values of the parame- 
ters, we trace U c in the parameter space. This is shown 
in Fig. 31 for Q/t = 0.5, 2, 4, 8 and various effective cou- 
plings. The MA results (filled symbols) are in good quan- 
titative agreement with the DMC results (empty sym- 
bols) for larger Q > 1 values (these are MA^ ) results for 
cutoff l c = 0. Using MA' 1 ' and/or increasing the cutoff 
changes the values of U c by less than 1% everywhere we 
checked). For the smaller frequencies such as ft = 0.5£, 
MA is known to become quantitatively less accurate at 
intermediary coupling a 15 ' 18 , and indeed, here we see a 
discrepancy with the DMC data even for the MA^ re- 
sults. To improve the quantitative agreement here, one 
should use a 2- or 3-site MA variational a ppr oximation 
for the phonon cloud, as discussed in Refs. [l7|. 

As expected, when A — > 0, U c goes towards the ex- 
pected critical value in the absence of electron-phonon 
coupling, of roughly 3.96£. As the effective coupling in- 
creases U c decreases, but the lines never intersect the 
y-axis: U c = is impossible, since the polaron cannot be 
trapped (localized) in a clean system as long as it has a 
finite effective mass, i.e. for any finite value of A. 

The decrease of U c with increasing A is expected, and 
is usually attributed to the fact that the effective polaron 
mass increases with A, and this makes it easier to trap 
near the impurityJ^ However, we claim that this is not 
the full story, and that the rcnormalization of the trap- 
ping potential also plays a non-trivial role. 

Consider, first, the Hamiltonian of Eq. (fT4")l , which de- 
scribes the impurity problem in the absence of electron 
phonon coupling. The binding energy of the impurity 
state (once formed) is a monotonic function of the only 
dimensionlcss parameter of this problem: -r- = / (y) 
for any U- > ^f- ~ 3.96. If one views the polaron as a 
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FIG. 6. (color online) Real part of the additional on-site 
MA' ' disorder potential vo(0, lu) over a wide energy range, 
for U = 2t, f2 = 2i, A = 0.8 and two values of 77. 



quasiparticle with an effective hopping t* which scatters 
on the same potential U as the bare particle (instanta- 
neous approximation), then the polaron binding energy 
should be |a = / (£.) for any f£ > £ « 3.96. In 
particular, this predicts U c /t = 3.96i*/t decreasing with 
increasing A, in qualitative agreement with Fig. 2) 

This hypothesis can be tested. The function f(x) is 
easy to calculate numerically, we can extract the bind- 
ing energy E* B for the trapped states by comparing their 
trapped energy to the GS energy of the polaron in the 
clean system, and the effective hopping t* ~ \/m* is 
directly linked to the effective polaron mass m* in the 
clean system.— We find that this scaling is not obeyed. 
Instead, one needs to also rescale the impurity potential, 

i.e. use -pr = f I ^7- J where U* =/= U. Of course, this 

scaling assumes that the scattering potential is local. As 
we discussed, this is not true although it is a good ap- 
proximation for medium and large couplings. 

In Fig. [5] we show the renormalized value U* /U ex- 
tracted this way, as a function of U/t above the cor- 
responding threshold values U c /t, for a medium and a 
large effective coupling A = 0.5, 1.5 and fl = 3i. Qual- 
itatively similar curves are found for other parameters. 
We see that U* — > U only when U — > oo and becomes 
the dominant energy scale (hence anything else is a small 
perturbation). For fixed U, fl, we find that U* increases 
with increasing coupling A - this is also expected, since 
the renormalization is directly caused by the electron- 
phonon coupling, see Fig. HJb). 

This renormalization is a direct illustration of the gen- 
eral result of Eq. (flU)) : the electron-phonon coupling 
changes not only the properties of the polaron (its effec- 
tive mass), but also the disorder potential it experiences. 

However, it is very wrong to expect that the poten- 
tial renormalization can always be described by a sim- 
ple rcscaling by some overall value. Indeed, Eq. (fl~3"|) 
shows that the renormalized potential is expected to be 
a function of energy, because of retardation effects. This 



function is not roughly constant, instead it has signifi- 
cant and very non-trivial dependence of u, as illustrated 
by plotting vo(0,uj) over a large energy range, in Fig. |6] 
Similar curves are found for other values of the parame- 
ters. Over a narrow range of energies around uj ~ — 7.5t 
where the bound state forms for these parameters (see 
Fig. [3]), uo(0,w) varies slowly and can be approximated 
as an overall negative constant. This explains why here 
we can approximate Cq(u>) = —U + Vo(0,u>) ~ —U*, with 
U* > U, as discussed above. At higher energies, how- 
ever, vo(0, ui) goes through singularities and changes sign 
from negative to positive and back. Although at first 
sight these singularities are surprising, they should be ex- 
pected based on Eq. (fTU]) . The self-energy of the Holstein 
polaron is known to have such singularities, especially at 
medium and higher couplings where an additional second 
bound state forms and the continuum above shows strong 
resonances spaced by fl. In particular, as A — > oo and 
the spectrum evolves towards the discrete Lang-Firsov 

2 

limit E n = — jy + nfl, the self-energy has a singularity 
at the top of each corresponding band. The renormal- 
ized potential of Eq. (fit))) is the difference between two 
such curves, displaced from each other. It is thus not a 
surprise that it has such nontrivial behavior. 

Physically, such strong retardation effects arc not sur- 
prising, either, since the additional potential vo(i,uj) de- 
scribes the scattering of the electron in the presence of 
the phonon cloud. The structure of the phonon cloud 
varies with energy, for instance one expects quite differ- 
ent clouds within the polaron band vs. at higher energies, 
in the continuum of incoherent states with finite lifetime. 
This suggests that the diagrams of Fig. HJb) that con- 
tribute most to the series change with energy, and so 
does the total result. As a final note, we mention that at 
these higher energies, MA' 1 ' should be used. It is well- 
known that MA' ' fails to describe properly the location 
of the polaron+one phonon continuum, since it does not 
include the needed variational states^ This problem is 
fixed at the MA' 1 ' and higher levels J£ 

To summarize, for this simple impurity problem the 
MA approximation is found to agree well with results 
from DMC in describing the trapping of the polaron. Al- 
though quantitatively not as accurate, besides efficiency 
its main advantage is that the analytic equations that de- 
scribe MA allow us to understand the relevant physics. 
In particular, we showed that coupling to bosons renor- 
malizes the disorder in a very non-trivial way. 



B. Impurity changing the electron-phonon coupling 

We now assume that the impurity does not change 
the on-site energy, but instead it modifies the value of 
electron-phonon coupling at its site: 
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FIG. 7. (color online) Phase diagram separating the regime 
where the polaron is delocalized (below the line) and trapped 
(above the line), as a function of the difference between the 
impurity and the bulk electron-phonon coupling, gd — g. Sym- 
bols show MA*- ' results, while the dashed black lines corre- 
spond to the instantaneous approximation. 



where gj = g + (gd — g)Sjo- Thus, gd and g are the 
electron-phonon couplings at the impurity site and in the 
bulk of the system, respectively. Since Vd = in this case, 
the non-interacting part of the Hamiltonian is fid = Ho- 
Thus, G d n {uj - nfl) -> Gf){u - nfl) = g (u - nfl) in 
the continued fractions, Eq. ^, whose dependance on 
the site index j is now through the coupling gj only. As 
a result, the effective disorder potential vq{j,u>) vanishes 
everywhere except at the impurity site, j ' = 0: 



vo(j,u) = A(w)5jo, 



(19) 



where A(w) = gaAi(Q, ui)— gAi(ui) , and A\(0,uj) is like in 
Eq. |9]) but with g —^ gd- This shows that even though 
€i = in this case, the inhomogeneity gives rise to an 
effective potential e*{ui) = 6ifiA(oj). This is now local 
because only when the phonon cloud is at the impurity 
site it can experience the different coupling. 

Equation pip can now be solved analytically to find: 



p(0,w) = Im 

7T \ 1 



ffo(w) 



A(u)g (u) 



(20) 



% = f + o y, h h + E wh fe + b \) 



(18) 



We now find the critical values gd — g when an impu- 
rity state emerges below the continuum, for given values 
of g,fl. The results are shown in Fig. [7] for gd > g, 
when the polaron formation energy at the impurity site, 
—g^/fl, is lower than the bulk value —g 2 /fl, and a bound 
state may be expected to form even within the instante- 
neous approximation. Symbols show MA' ' results, while 
the dashed lines are for the instantaneous approximation. 
The two agree quantitatively only in the limit A — > 0. 
This proves that the additional renormalization included 
in MA is significant for this type of impurity, as well. We 
note that all critical lines intercept the rr-axis at a finite 
value, i.e. for any value of fl and g = 0, there is a critical 
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FIG. 8. For large A, in the clean system (dashed red line) the 
first polaron band is separated by an energy gap from the next 
features in the spectrum (here, the band associated with the 
second bound state). For gd < g, an "anti-bound" impurity 
state is pushed inside this gap (black full line). Parameters 
are fi = t = 1, A = 1.2 and r\ = 10" 3 . 



interpret this as being an "anti-bound" polaron state. A 
similar state is also expected to appear for a repulsive 
on-site U < potential. 



C. Isotope impurity 

The last case we consider in detail is an isotope im- 
purity. Because of its different mass Md ^ M, both its 

phonon frequency fid ~ M. 2 , and its electron-phonon 

_ i 
coupling gd ex 1/ \JMdfld ~ M d 4 , are changed. Inter- 
estingly, both the effective coupling Ad = g d /(Qtfld) = 
g 2 /(6tft) = A and the polaron formation energy 
— g^/fid = —g 2 /Vt show no isotope effect^ As a result, 
within the instantaneous approximation of Eq. (|12[) one 
would predict that the isotope is "invisible" and the po- 
laron spectrum is basically unaffected by its presence. 

We consider a single isotope impurity located at the 
origin: 



finite value gd above which an impurity state forms. For 
example, for D, = 2>t this critical value is gd ~ 3.9i. Its 
value increases with increasing $7, as expected. 

Unlike for an impurity which changes the on-site po- 
tential, and which can bind at most one impurity state, 
impurities which change the electron-phonon coupling 
can bind multiple impurity levels. As gd increases and 
the energy of the impurity level moves towards lower 
energies, additional bound states, spaced by roughly fi, 
emerge whenever the distance between the last one and 
the bulk polaron band is of order $7. 

The origin of this sequence of bound states is straight- 
forward to understand in the limit gd ~> g,t, where, the 
Hamiltonian is, to zero order: 

H « g d clc {b\ + b ) + £!6 6 , 

with c cq ~ 1 because the weight of the bound state 
is concentrated at the impurity site. This Hamilto- 
nian can be exactly diagonalizcd with the Lang-Firsov 
transformation^ and predicts a series of equally spaced 
eigenenergies nil — g^/fi. For finite t, g, all states which 
lie below the bulk polaron continuum become impurity 
states, and basically describe excited bound states with 
additional phonons at the impurity site. 

So far we have considered gd > g, where a ground-state 
impurity level can emerge. It is important to note that 
discrete peaks can also appear for gd < <?, although not 
at low energies. This happens when A is sufficiently large 
that there is a gap between the bulk polaron band and 
higher features in the spectrum, such as the polaron+one 
phonon continuum, or the band associated with the sec- 
ond bound state, once it forms^ A typical example is 
shown in Fig. ___, where in the presence of an impurity 
with a weaker coupling gd < 9 (full line) , a discrete state 
appears above the polaron band. Since most of its weight 
is removed from the bulk polaron band (dashed line), we 



*■ = t + E whih + b]) + _T ntfbj, 



(21) 



where fL 



fi 



"3 — "T (fid 

chosen such that Aw 



9 + (9d - g)Sjo are 



fi)<5j0 and g 3 - 
A. 

Just like in the previous section, because there is no on- 
site disorder e, = 0, we have G^{uj) = Gj (w) and the 
effective disorder potential is again local, i.e. it vanishes 
everywhere but at the impurity site. As a result, the 
LDOS at the impurity site is given by Eq. (J____|), except 
that here: 



A(w) = E d (w) 



J MA< 



o)M, 



(22) 



where S^(w) has the same functional form like S MA <o) (w), 
but with g — > gd, fi — >• fid- 

From investigations of this LDOS at different param- 
eters we find that there exists a threshold value of the 
effective coupling, A*, below which low-energy bound 
states do not form irrespective of how small or large 
Md/M is. In other words, for A < A*, the behavior agrees 
with the prediction of the instantaneous approximation. 

We can estimate a bound on A* as follows. Consider 
the case of a very light isotope, so that fid,5d S> fi,ff,t- 
In this limit, £<i(u;) — > —g^/Cld = — 6tA. The bound state 
appears when the LDOS is singular because its denomi- 
nator vanishes: 

1-A(uj)g (u>) = -> S MA( o, (w)+l/<7 (w) - -6tA, (23) 



after using Eq. ([22]) and Sd(tf) ~ —6tX. 

Consider now the limiting case when a bound state 
emerges just below the bulk polaron ground state, i.e. 
Eq. (f2"3"]l has a solution at u> < e p ° l . In the clean system, 
the polaron ground-state energy e^f is the lowest pole of 
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p pol 



0,w) 
-6t+ 



£k=0 
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£ MA (o)M]~ 



so it satisfies: 



S M a (o) (e p ° l ) . Using this in Eq. (f2"3"| suggests 
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FIG. 9. (color online) Critical effective coupling A* above 
which an impurity state may appear for a sufficiently light 
isotope. Below this line, polarons cannot be bound near iso- 
topes. Symbols shows MA' ' results. The dashed line is the 
analytic low-bound for A* discussed in the text. 



that a solution can exist if A > A* , where 



,pol 



(i/ 



- 1 - 



1 



6i 3o (-6t) : 



(24) 



with ff (-6t) ~ -1/3.96*. Since s p g f -> -6i as A ->■ 
0, we expect that A* — > 0.66 in this limit, and that it 
increases as e^° l becomes more negative, for example with 
increasing A. These considerations are confirmed by the 
data shown in Fig. [5] Here, the symbols show values of 
A* found numerically with MA' ) , and the dashed line is 
the lower bound of 0.66, discussed above. 

For A > A*, bound impurity states can appear near 
isotopes if gd and fid = g^/(6t\) are sufficiently large. In 
Fig. [TU] we show critical lines for two cases, fl = 4t, 8t. 
The symbols show the MA' ) results, which converge to- 
wards their corresponding A* values as fid — > oo, as ex- 
pected. Of course, the largest values considered for fid 
are unphysical; we use them only to illustrate the con- 
vergence towards A*. 

The existence of a region of the parameter space where 
bound polaron states appear near an isotope is in direct 
contradiction of the instantaneous approximation, and 
again illustrates the importance of the rcnormalizcd dis- 
order potential A(w) which makes their trapping possi- 
ble. In this context, it is worth mentioning that there 
is clear evidence for electronic states bound near isotope 
O 16 defects in Cu02 planes^ although the precise na- 
ture of these states has not been clarified and the mea- 
surements are certainly not in the extremely undcrdoped 
regime where our single-polaron results are valid. 

Interestingly, when such bound states form near an 
isotope, the spectrum is different than that for the other 
two types of impurities. As shown in Fig. II H bound 
states now appear simultaneously both below and above 
the bulk polaron continuum, not just below it. This pro- 
vides a possible "fingerprint" for polarons trapped near 
isotopes. Finally, we note that even when no low-energy 
impurity state is observed, it is again possible to have 



higher energy bound states inside the gaps opening be- 
tween various features in the bulk polaron spectrum. 

To summarize, in the presence of isotope defects, po- 
larons in the weakly coupled regime A < A* always re- 
main delocalized. Only for A > A* it is possible to trap 
polarons near an isotope. This makes this case quite dis- 
tinct from the other two cases, where bound states exist 
for any A if the impurity is strong enough. 



IV. SUMMARY AND CONCLUSIONS 

We studied the threshold for the emergence of polaron 
bound states near various types of single impurities, using 
the momentum-average approximation. Elcctron-phonon 
coupling was shown to strongly renormalize the impurity 
potential in a nontrivial way which includes strong retar- 
dation effects. This is a feature that is completely absent 
in the instantaneous approximation, which is the only 
other available "simple" description of this problem. 

We considered the simplest models of impurities that 
change the strength of the on-site energy, the local 
electron-phonon coupling, or arc isotope substitutions 
that modify both the coupling and the phonon energy. 
We calculated the polaron binding phase diagrams for 
each case. The first case had been considered previously 
by numerical methodS ) 12 i 13 and our results are in good 
quantitative agreement with their predictions. To our 
knowledge, the other two cases have not been investi- 
gated before. We showed that in the first two cases 
bound states always exist for a sufficiently strong im- 
purity, however the polaron remains delocalized for the 
case of isotope substitution of arbitrary strength if the 
effective coupling is weaker than a threshold value, A*. 
Differences in the LDOS at the impurity site have also 
been found, such as the possibility to bind multiple states 
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FIG. 10. (color online) Phase diagram separating the regime 
where the polaron is delocalized (below the line) and trapped 
(above the line), as a function of the difference between the 
fid - fi, on a logarithmic scale. Symbols show MA' ' results 
for fl = 4t, 8t. As fid —¥ oo, these critical lines converge 
towards their corresponding A* (dashed lines), below which 
polaron states are always delocalized. 
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FIG. 11. (color online) LDOS at the impurity site near an 
isotope with fid = fi + It (full line). Two discrete states, one 
above and one below the bulk polaron band, are seen. For 
comparison, the DOS in the clean system (multiplied by 10) 
is shown as a dashed line. Parameters are fi = 4t, X — 2.5 
and r\ = 1CT 2 . 



near an impurity that changes the coupling, or the un- 
usual fingerprint of discrete states both below and above 
the bulk polaron continuum, for an isotope bound state. 

Of course, a realistic description of an impurity in 
a real system may combine several of these inhomo- 
geneities, and even the form of the electron-phonon cou- 
pling could be affected. MA gives an efficient yet quite 
accurate way to deal with such cases, and can be easily 
generalized to other types of couplings where MA has 
been used succesfully to describe bulk properties. 

Whereas we expect the single impurity results to re- 
main valid for a system with multiple impurities if the 
mean free path is long and the polaron interacts with one 
impurity at a time, in the presence of significant disorder, 
when multiple scattering processes become important, 
the polarons can undergo Anderson localization. This 
limit has been addressed within a generalized DMFT^. 
however we believe that our simpler formulation might 
provided additional insight and uncover previously un- 
explored aspects of Anderson localization for polarons. 
Such work is currently in progress. 
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Appendix A: IMA (1) 

At the IMA' 1 ) level, one also allows processes in 
which one phonon is away from the phonon cloud. 



(0\ciG(u)cJbj n ~ b]\0) with j ^ I. In terms of these, Eq. 
([3]) can be written as 

G«(«) = g£-(w) + £>s 1 (»,M;w)G«M- ( A1 ) 



Once again we apply the Dyson identity to Si 

Si (i, I, j; w) = g 3 G% (w - fy )G ij (w) (A2) 

+ ^2gmG d ml {uj - flj)S 2 (i,m,j;u). 



This exact equation relates Si to the propagators 52- 
We can similarly find the equation of motion of all higher, 
n > 2, S n (i,l,j;u>) with I ^ j and I = j, separately. For 
I =/= j we have 

S n (i, l,j\u)= giGf^u - (n - 1)0, - fy) (A3) 

x[(n- l)S n -i(i,l,j;w) + S n+ i(i,l,j;u)}, 

where we now ignore contributions from terms with a sec- 
ond phonon away from the polaron cloud, as they are ex- 
ponentially smaller than those we kept. This admits the 
solution S n (i,l,j',oS) = B n (l,j;u)S n -i(i,l,j',w), where 



Bn(l,j;u) 



(n - ^giGftJu - {n - l)Qi - Qj) 
1 - giGf^ - (n - l)fi, - Q,)B n+1 (l, j;u) 
A n _i(l,u;-Q j ). (A4) 



For I = j and n > 2, S n (i,l,l;u) = F^ (w) and we 
get the same solution as in MA' ', i.e. S n (i,l,l;oj) — 
A n (l,oj)S n -i(i,l,l',u)- The relations between 52 and Si 
are used in Eq. (|A2j) to turn it into an equation between 
Si(uj) and Gij(to) only 

Si(i,l,j;u) = gjGi{u - Q,j)G io (w) 



+ 9jGji (w - toj)A 2 (j, w)Si (i, j,j;u) 
+ J2 gmG^tiu - nj)Ai(m,u - ttj)Si (i,m,j;w). 



E 



Together with Eq. (|A1|) . this can be solved to find 
Gij(ui). However, it is again convenient to explicitly ex- 
tract the "average" contributions, to make these equa- 
tions more efficient. 

We therefore remove the homogeneous part from Eq. 
[2| and include it into a renormalized energy: 



12 



Si(i,l,j;u) = g 3 G%{uj 3 )[G l3 (uj) + (A 2 (j,u) - A x {j,u - ^))Sx{i,j,r,u)} 

+ Y. 9rnGti(uj)[Ai(m,u - fi,-) - Ai(w - n)]Si(*,m,j;w). (A5) 

where Wj — w — ilj — gAx(u — 51), and A\(u — ft) is given by Eq. ([9]) for the "average" clean system. 

The sum on the rhs of Eq. (|A5j) again convergences for a very small cutoff, only sites m very close to j need to be 
included. Its general solution is of the form: 



where 



Si(i,l,j;uj) = Xji(u)[Gij(u) + (A 2 (j,u) - A x {j,u- flj))Si(i,j,j;u)], 






(A6) 



In fact, using x 3 i(uo) = gjGJ^uij) is already a very good efficient to solve by subtracting the MA^ self-energy 
approximation, since the terms in the sum are expo- and including it into the energy argument 
nentially small - but one can go beyond this. Once 
Xjjitu) is known, from Eq. (|A6|) we find Si(i,j,j,u>) — 
Aj(u))Gij(ui), where 



. / \ x 33 (uo) 

AU) ~ 1 - x 33 (uj) {Mj, uj) - Mj, uj - Qj)) ■ 

This can now be used in Eq. (jAll) to turn it into an 
equation for Gi 3 {w) only: 

G«(w) = Gfj(u) +Y,g l A i (u)G a (w)Gt j (u). (A7) 



Gij(u) = GijW + J^Gui^Mh^GUu), (A8) 
i 

in which uj = uj — S MA <i) (uj) and v\(l,oS) = giAi(u>) — 
E MA (i)(w). Here, S MA (i)(w) is the value of giAi(uj) in 
the clean, "average" system: 



J MAW 



(«) 



3 2 3o(w) 



1 - .g 5o (w) (A 2 (w) - j4x(w - fi)) 



where now <D = w — Jl — (/j4i(w — SI). 18 This completes 
the calculation of Green's function within inhomogeneous 
As we did in Eq. (fTT|) for MA' ), this can be made MA^ 1 ' approximation. 
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